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Depth averaged conservation equations are written for granular surface flows. Their application 
to the study of steady surface flows in a rotating drum allows to find experimentally the constitutive 
relations needed to close these equations from measurements of the velocity profile in the flowing 
layer at the center of the drum and from the flowing layer thickness and the static/A owing boundary 
profiles. The velocity varies linearly with depth, with a gradient independent of both the flowing 
layer thickness and the static/flowing boundary local slope. The first two closure relations relating 
the flow rate and the momentum flux to the flowing layer thickness and the slope are then deduced. 
t-H ' Measurements of the profile of the flowing layer thickness and the static/flowing boundary in the 

whole drum explicitly give the last relation concerning the force acting on the flowing layer. Finally, 
. these closure relations are compared to existing continuous models of surface flows. 
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! I. INTRODUCTION 

T— I ! 

^ • Granular media present numbers of interesting and unusual properties |l| : the intrinsic dissipative nature of the 
' interactions between the constituent macroscopic particles sets granular matter apart from conventional gases, liquids 
or solids. One of the most interesting phenomenon in granular systems is the transition from a static equilibrium to 
a granular flow: contrary to ordinary fluids, they can remain static even with an inclined free surface. But when the 
angle of the surface exceeds some threshold value 6^, the pile can not sustain the steep surface and starts to flow until 
its angle relaxes under a given angle of repose 8 r . The motion has the particularity to be a surface flow: most of the 
pile remains static. The condition governing the transition between the solid state and the liquid state as well as the 
internal equations of the flowing layer are still today under debate 0. This makes the description of an avalanche 
rather difficult. 

People have benefited from the particularity of avalanches to be surface flows to apply to them depth averaged 
conservation laws, the Saint Venant equations ||], and to use the approximations developed in hydrodynamics for 
a thin fluid film flowing down an inclined plane. Savage and Hutter JiJ], and more recently Pouliquen [5 applied 
these equations to describe the motion of granular material down a rough inclined plane. Khakhar et al. [3|, and T. 
Elperin and A. Vikhansky |t| have studied steady flows in a rotating drum using the same theoretical framework. 
They succeed to report on the different free-surface shapes successively observed when the rotation speed increases. 
Finally, Douady et al. ^ and very recently Khakhar et al. Q have used these equations to describe avalanches at the 
surface of a heap. Another model, referred to as the BCRE model [|l0|, has also been proposed a few years ago: it 
deals with thin surface granular flows including avalanches and succeeds for example to report on the existence of the 
two angles 9d and 9 r . This model has then been extended to describe thick sur face flows 0Q. 

All these models are based on the assumption that a strict separation between flowing grains and static grains 
can be made. They end in a system of two coupled equations with the same form, differing only through their 
closure relations. Those are chosen from theoretical considerations, linking for example the coupling term to the 
effective force applying on the flowing phase (for Saint- Venant approach) or to the microscopic effects of the flowing 
grains collisions on the static bed (for BCRE approach). An original approach, proposed by Aranson and Tsimring, 
should also be mentioned JT3]. Assuming a continuous transition between flowing grains and static grains, they have 
developped a theory based on hydrodynamic equations coupled with an order parameter equation describing the 
solid/flowing transition. This theory accounts for most of the behaviours recently observed in experiments on rough 
inclined plane (T^| . 

Thin rotating drums, i.e. rotating drums whose gap is small compared to the diameter, have been widely used 
to investigate properties of granular materials: Evesque and Rajchenbach jl5j and Jaeger et al. (TsJ have studied 
avalanche size distribution. Rajchenbach has used it in an attempt to find constitutive mechanical laws in 2D beads 
packing |p^| . Nakagawa et al. |]l8f have used magnetic resonance imaging to measure density and velocity field in 
a thick rotating drum. We propose here to study steady surface flows in a rotating drum for a quasi-2D packing 
of monodisperse beads in order to find experimentally the closure relations needed to complete the depth averaged 
hydrodynamic equations. 
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In section |l[ basic conservation laws are presented. It is shown how the closure relations can be deduced from 
the experimental velocity profile inside the flowing layer at the center of the drum and from the flowing layer profile. 
The experimental set-up is then described (section HI). Beads are studied on the microscopic scale in section IV: 
velocity and volume fraction profiles are measured on elementary slices via a fast camera. In section 0, the flowing 
layer thickness and the angle between the interface and the horizontal are measured in the whole drum. The closure 
relations are deduced. Finally, they are compared in section [v| to the ones assumed by the different models proposed 
in the litterature. 



II. THEORETICAL FRAME: SAINT- VENANT DESCRIPTION 

Let us consider the situation depicted fig. The approach is limited to the 2D situation. The material can then 
be divided into two parts: a thin cascading layer flows on a static bed. The boundary between these two phases is 



supposed to be sharp. This assumption will be discussed and justified in section |VJ. As the flowing layer is assumed 
to be thin, the velocity of grains is essentially parallel to the static/flowing boundary. The most natural frame is 
consequently the curvilinear frame (et, et) where (resp. et) is locally tangent (resp. normal) to this interface. 
The z origin is set at the free surface. Calling (v x , v z ) the velocity at the (x, z) coordinate, the mass and momentum 
conservation equations can be written 



— (< p > R) + —(R < pv x >) + (pw,)U=_ fl = (1) 
d d 

— {R < pv x >) + —{R < pv x 2 >) + {pv x v z )\ z= - R = F (2) 

where R is the flowing layer thickness, p the solid fraction and F the x-component of the force acting on the volume 
of the flowing layer bounded by the two surfaces £(x) and E(x + dx). The quantities (R < pv x >) and (R < py x 2 >) 
are respectively the flow rate and the flux of x-componcnt of momentum through the surface S(cc) (see fig. |f|) while 
(pv z )\z=-R and (pv x v z )\ z =-R represent respectively the flow rate and the flux of x-component of momentum passing 
through the static/flowing boundary. < A > denotes the average of the quantity A across the flowing layer : 

1 r° 

< A >= — / Adz (3) 

R J-R 

In the following, the volume fraction p is supposed to be constant (this assumption will be discussed and justified 
in section [V) and eq. (QJ^) are simplified accordingly. The expressions of < v x >, < v x 2 > and F as a fonction 



of the flow macroscopic parameters are required to close the two equations (|1|,[2|). Translating the packing upwards 
should not change its dynamical equations. Consequently, the three closure relations can be written as < v x > (R,9), 
< Vx 2 > {R,&) and F(R,9) where 8(x,t) is the angle made by the solid/flowing interface with the horizontal. Let 
us apply these two conservation equations to the description of steady surface flows in a rotating drum (see fig. ^). 
In this case, time derivatives terms in eq. vanish and the velocity v\ z =-r of the static grains at the bed/layer 

interface can be written: 

v\z=-R = ttrn (4) 

where is the rotation speed of the drum, r the distance between the center of th e drum O and the point M of the 
interface where v\ z= -r is calculated, and ft the unity vector perpendicular to oM. Then, eq. (0,0) become: 



d (R < v x >) = -n^Xcose + HsmeJ (5) 



dx 



-^(R < v 2 >) + n 2 ^X cos 9 + H sin 6>) (x sin 6 - H cos 0) .H (x cos 9 + H sin 6j = ^ 



where (X(x),H(x)) is the coordinate of M in the gravity frame whose origin coincides with the center of the drum 
(see fig. 0) and H the Heaviside function. 

The aim of the paper is first to investigate the beads behaviour on the microscopic scale: the profiles of the velocity 
v x (z) and volume fraction p{z) are measured on an elementary slice. The evolution of these profiles with R and 9 
allows one to deduce the first two closure relations < v x > (R,9) and < v x 2 > (R,9). Then, profiles R(x) and 9(x) 
are measured on the macroscopic scale, in the whole drum. These profiles confirm the form of the first two closure 
relations, extend their range of validity and allow one to find the last closure relation F(R, 9). 
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III. EXPERIMENTAL SETUP 



The experimental set-up is illustrated in fig. g. It consists mainly in a duralumin rotating drum of diameter 
Do = 45 cm and of gap e = 7 mm. The drum is half filled by a packing of steel beads of diameter d = 3 ± 0.05 mm. 
The drum can then be considered as a Hele Shaw cell (Do > e k 2d) where the packing has a 3D microscopic structure 
(especially a disordered structure contrary to 2D packing whose stable structure is crystalline) while keeping a 2D 
macroscopic geometry. The volume fraction measured just after the filling is 0.57 ± 0.02. The gap e has been chosen 
slightly larger than 2d to prevent any jamming in beads flow. The use of steel beads of millimetric size allows to have 
a good control of geometrical and mechanical properties and to limit capillary and electrostatic effects at the expense 
of a larger scale experiment. 

The rotation speed f2 can be varied from 0.001 r.p.m. up to 30 r.p.m. For small Q - typically smaller than 0.1 r.p.m. -, 
the mean angle of the free surface 9 oscillates between the angle of flow 9 c i = 27.5 ± 0.7° and the angle of repose 
r = 23.5 ± 0.6°. This intermittent behaviour is not studied in this paper. For f2 > 0.1 r.p.m., surface flow becomes 
steady. The beads are lighted via a continuous halogen lamp and regimes obtained for Q varying from 1 r.p.m. up 
to 6 r.p.m. are recorded via a fast camera. On this range of rotation speed, the width of the flowing layer R is well 
defined and inertial effects are negligible (the Froude number Fr = n 2 D /2g < 0.01 where g is the gravity constant). 



IV. MICROSCOPIC SCALE: VELOCITY AND VOLUME FRACTION PROFILE ON AN ELEMENTARY 

SLICE 

Sequences of frames are recorded via a fast camera at a sampling rate f s = 1 kHz. The shutter speed is of 0.1 ms. 
The recorded region is located at the point where (vt\ z =-R- ~ 7 ft) = (roughly corresponding to the center of the drum). 
In this region, 9, R and {v x (z)} are invariant under a small translation along ~e* x . The recorded region can then be 
considered as a statisfactory elementary slice. R and 9 in this area are controlled by changing Q, but can not been 
modified independently. The recorded region size in pixel is 480 x 234, one pixel corresponding to 0.227 mm. The 
digital images are processed to obtain the position of the center of mass of the beads seen through the glass porthole 
(see fig. ||). The image of a single bead is made up of about 23 ± 3 pixels depending on the distance of the bead to the 
porthole. The error on the determination of beads location is thus about 50 /xm. One can then extract both volume 
fraction and velocity profiles. 

Volume fraction profile can not be measured directly in our experiment since all beads are not seen through the 
porthole. One has thus calculated from numerical simulated packing the relationship between the real volume fraction 
and the volume fraction estimated from the beads seen through the porthole (see appendix A). The experimental 
volume fraction profiles can then be deduced. The errors made on these profiles is mainly dominated by the statistical 
processing of the data: profiles have been averaged over 1000 frames to limit this error. A typical volume fraction 
profile is represented fig. |ja: p varies typically from 0.3 to 0.4 in the flowing layer, except at the free surface in a small 
region (about three beads diameters) where p drops quickly down to zero. This last region is assumed to be sharp 
and corresponds to the free surface boundary. The assumption p ~ constant can then be questioned: relaxing this 
assumption amounts to replace < v x > and < v x 2 > respectively by < pv x > / < p > and < p(z)v x 2 > / < p >. As 
p variations are negligible compared to v x variations, p can be assumed to be constant. 

Frames processing allows to evaluate beads velocity (see appendix B). A typical velocity profile is represented fig. [|b. 
The velocity profile is linear in the whole flowing layer. The inset, representing the variation of the velocity vs. depth 
in semilogarithmic scale, emphasizes the creep motion of the static phase. In this phase, v x decreases exponentially 
with the depth with a characteristic decay length A creep = 7.6 mm ~ 2.5d in agreement with p9| . The creep motion is 
thus localised in a narrow layer at the solid/liquid interface whose thickness can be considered as zero at macroscopic 
scale. The creep motion is not studied in this paper. 

The rotation speed f2 has been varied from 1 r.p.m. to 6 r.p.m. in order to study variations of the velocity profile 
with R and 0. The corresponding velocity profiles have been represented fig. |B|a. In fig. ||b, all profiles have been 
translated along the z-axis to make the static/flowing interface coincide. The velocity profile remains linear in the 
whole flowing layer and the velocity gradient is independent of both R and 0. 

Velocity profiles measured at five different points equally distributed along the free surface for O = 6 r.p.m. are 
represented in fig. ^. At these locations, d x R and d x 9 are no more equal to zero. The thickness of the recorded region 
has been decreased to minimize these drifts. {v x (z)} does not depend on the location in the drum, which means that 
{v x (z)} does not depend on d x R and d x 9. The velocity profile can be then written inside the flowing layer as: 

v x (z) =j(z + R) (7) 
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with 7 = 34± 0.5 s _1 , or 7 ~ 0.6-\/<7/d. Such order of magnitude can be easily understood p7[ : when a bead collides 
with an underlying grain, all the kinetic energy is absorbed by multiple collisions. Then, the balance between potential 
and kinetic energy on a mean free path of order of d leads to a lim iting velocity between two adjacent beads layers 
equal to y/gd, and consequently to a velocity gradient of order of \ fgjd. 

The knowledge of the profile {v x (z)}(R, 6) enables us to calculate the first two closure relations: 

<v x > (R, 6) = ^7-R (8) 



< v 2 > (R,9) = p 2 R 2 (9) 

This scaling law is different from the one observed in granular flows down rough inclined plane for 2D packing ]2G[ ] 
and for 3D packing [pT| : I n these experiments, < v x > was found to scale as i? 3 / 2 . A linear profile with a velocity 
gradient independent of R has also been observed by Rajchenbach at the surface of a 2D packing in a rotating 
drum ]Ti| . This tends to confirm that boundaries conditions can strongly affect the velocity profile and questions the 
ability of a classical constitutive law relating the shear stress to the strain rate PJl7||. 

Let us note that velocity gradient observed by Rajchenbach in 2D packing scales as Vsin0 contrary to ours. However, 
in our series of experiments, 9 varies only from 31° up to 48°, which corresponds to variations of V sin 9 from 0.72 to 
0.86. Consequently, the difference between 7 behaving as Vsin^ and 7 constant does not differ much from the error 
bars and it is not possible to conclude at this point which of the two scaling is the right one. This will be tested on a 
more significant range of 9, from 0° up to 48°, in the next section. 



V. MACROSCOPIC SCALE: PROFILES OF THE FLOWING PHASE HEIGHT AND OF THE 

STATIC/FLOWING BOUNDARY 

The determination of the last closure relation F(R, 9) requires the determination of R and 9 in the whole drum. 
The experimental protocol is the following: sequences of 512 frames arc recorded via a fast camera at a sampling rate 
f s = 250 Hz with a pixel resolution of 480 x 400. In order to keep a sufficient resolution on the profiles, the rotating 
drum is visualized by halves, a pixel corresponding to 0.468 mm. Image processing allows to isolate the flowing layer 
and both the interface profile H(X) and the free surface profile S(X), expressed in the gravity frame (ex,ez), can 
be extracted (cf fig. [?]). The fluctuations on each point of these profiles are around 1.5 mm. The profiles are then 
averaged over 20 frames. Consequently, the resulting errors on both H and S are around 0.3 mm. The curvilinear 
coordinate x(X), the local angle 9(X), and the flowing layer thickness R(X) are then deduced using the procedure 
given in appendix C. Resulting errors on R is about 1 mm (see fig. ||b). On the upper part of the profile, the local 
slope is constant on a range whose size decreases when f2 is increased (see fig. ||a) . One can study on this range Cq 
the variation of the different variables with R for a fixed angle 9q depending on the rotation speed fi. 

Let us first confirm the first closure relation (eq. ||) established at the end of the last section: the drift velocity 
< v x > can be calculated from the experimental x(X), H{X) and R(X) profiles by integrating numerically eq. [jj. 
Fig. ||a shows the variation of < v x > vs. 57-R with 7 = 34 s _1 as suggested in the previous section for different 6q. 
This confirms the assumption of a linear profile {v x (z)} with a velocity gradient independent of R. 

Let us answer now to the question about the dependence of 7 with 9. To test the different assumptions on this 
scaling, namely if 7 scales as Vsin# or is independent of 9, profiles H{X) and R{X) have to be taken into account in 
their whole: Typically for £1 — 6 r.p.m., 9(X) varies from 0° to 48°. The procedure will be the following: from the 
profile H exp (X) obtained experimentally one deduces X(x), H(x), 9(x) and, by integrating the eq. the flow rate 

{R < v x >}(x). The theoretical profiles for the flowing thickness R^(x) and R^(x) are then calculated assuming 
respectively the relations: 

(a) R < v x >= pR 2 corresponding to a velocity gradient independent of 9 

(b) R < v x >= pR 2 Vsm9 corresponding to a velocity gradient 7 scaling as v / sin# as observed by Rajchenbach (Tt]] 

These two profiles are then compared to the experimental profile R exp (x) (see fig. ^b). The fit parameter is p. The 
scenario (a) leads to better fits, with a fit parameter p — 14 s" 1 about 20% smaller than 7/2. This difference may 
come from the boundaries: the velocity gradient 7 has been determined from the velocity of only the beads observed 
through the porthole while p is related to the mean velocity of all the beads. 
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At this point, one can evaluate the force F. Numerical calculation of the second left-hand term of eq. ^ shows that 
it is negligible compared to the first left-hand term and can thus be neglected. From eq. the simplified eq. |6| and 
the closure relations || and |[ one finds: 

— = h sine + x cos e) (10) 



p 

The last closure relation can then be determined. From the experimental X(x), H(x) and R{x) profiles, one can 



deduce F(x) using eq. 10, Fig. |lO| a presents the variation of F/fryR vs. R for different 9. F/pjR seems to decrease 
as a parabola with R: 

F x (6) - F 2 (9)R 2 (11) 



piR 

A minimal form of the force F should include: 

• the weight Fyy — pgRsmO 

• a pressure term Fp = —kpg^ cos 9 Since R is negligible compared to D , this term can be neglected. 

• the friction on the static phase: Fs.fr = —ppgRcos9. A priori, p depends on R and 9. Experiments shows 
however that p reach quickly a constant value close to tan 9 r independent of R and 9 as soon as R is larger to 
few beads diameters p^ |. In our case, p, can be considered as constant. 

• finally, the friction on the lateral boundaries Fs.fr = Jq 'cs./rdz. 

This last contribution is the one responsible for the decrease of F/pjR with R as a parabola. It can be obtained 
directly from the value of F and the first two contributions Fw and Fs.fr- Fig [Io| b represents the variation of Fs.fr / p 
vs. R for different 9n. Its variation with On appears to be weak. The form of Fs.fr can be understood through the 
following process: Fs.fr = / &B.frdz where as.fr(z) is the stress exerted by the boundary on the layer between 
z and z + dz. To model Of r .s, let us apply the following process: Ofr.B = v8p where v is the frequency of bead- 
boundary collision and bp the x-momentum lost by each collision. As p is constant, the mean free path is constant in 
the layer and v oc pv x (z). Each collision with the boundary is inelastic. Consequently, dp cx v x . Finally one can write 
cfr.B — f]P v x 2 with r] constant and, after integration over the width of the flowing layer: 

. 2 

Ffr.B = -P^-R 3 (12) 

This closure relation indeed allows one to reproduce the experimental Fs.f r (R) profiles obtained for different 9q 
(cf. fig |lO[b). In our experiment, 77 ~ 10 m — 1 . The order of magnitude of 77 can be understood in the following way: 
the momentum lost by each collision is of order (1 — e) where e is the restitution coefficient. For our beads, e ~ 0.95. 
The mean free path is of the same order of the bead diameter d. Consequently, one expects 77 ~ (1 — e)/d ~ 10 m . 

The last closure relation that follows from the above discussion thus reads: 

F = pgR sin 9 - ppgR cos 9 ~ ^r^i? 3 (13) 



VI. DISCUSSION 



Using the closing relations eq. ||, eq. ^ and eq. 13, the conservation equations (QJ^) become 



d_R + v d_R =r 

dt dx 



dH 
~dt 



= -r 



(14) 
(15) 



where: 
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V = 2<v x >= 7-ff 



(16) 



and: 



v(sin0- ucos(9) - -n-viT (17) 
/57-K 7 ' 3 

This system of equations contains three parameters: the velocity gradient 7, the friction coefficient /j and a third 
parameter 77 characterizing the friction on lateral boundaries. These three parameters depend on microscopic features 
of the beads such as their diameter, their coefficient of restitution and their roughness and of the gap of the drum. It 
can be compared to the different existent continuous models proposed in |§-l3]. 

These two equations are very close to the one proposed by Douady et al. \ W. They have supposed a linear velocity 
profile with a constant velocity gradient as observed in our experiment. Moreover, for large R, variation of the force 
F with R and 9 is qualitatively similar to ours, but they attribute this variation to the friction of the flowing layer to 
static grains instead of invoking friction effects on the lateral boundaries as we have done. 

Khakhar et al. Q have proposed a derivation of depth averaged conservation equations very close to that of [||. 
They also have assumed a linear velocity, but with a gradient dependent of 9. As for the closing relations proposed 
by g, variation of F with R is attributed to a dependency of the shear stress at the solid/flowing interface with R 
rather than to boundary effects. 

Another model has also been proposed a few years ago by Bouchaud et al |l(J. The original model, devoted to 
describe avalanche triggering, deals with thin moving layers of R < d and can not be compared to ours. But this 
model has later been extended to describe thick surface flows by Boutreux et al. [jll| . These models assume a constant 
advection velocity V, small variations of 9 around 9 r , and relate the coupling term T to the effect of collisions between 
the moving grains on the static bed. For thin flowing layers [HJ , all the moving grains are assumed to be in contact 
with the static bed and to have the same probability to interact with it: T = aR(9 — 9 r ) whith a constant. For thick 
granular flows [jy], only grains belonging to the first layers can interact with the static bed. The influence of the 
others is shielded by these first layers. T is then more naturally assumed to be independent of R: T = a(9 — 9 r ) 
with a constant. This last form is compatible with our form of T expanded to first order in 9 — 9 r . The constant 
a can be evaluated ||J|: a — gcos9 r /j. Thanks to their simplicity, these two last model allow analytical exact 
solutions p3|-p5| . But the assumption of a constant velocity V prevents them for describing quantitatively avalanches 
front p^j or avalanches amplitudes fl2]j . 



VII. CONCLUSION 



Depth averaged conservation equations have been written for granular surface flows. Their applications to the study 
of steady surface flows in a rotating drum allow one to find experimentally the closure relations of these equations from 
measurements of the velocity profile in the flowing layer at the center of the drum and from the flowing layer thickness 
R and the interface slope 9 profiles. These measurements have been performed for a quasi 2D packing of steel beads. 
The velocity profile has been found linear with a gradient 7 independent of both R and 9. The closure relations 
have been deduced: the mean drift velocity < v x > and the mean kinetic energy < v x > are given respectively by 
< Vx >— \iR and R < v x 2 >= ^j 2 R 2 . The last relation relating the force F acting on the flowing layer to R and 9 
has also been explicitly given and reveals the importance of lateral boundaries. Two parameters, namely the friction 
coefficient /1 of the flowing layer on the static bed (close to tan0 r ) and a coefficient 77 characterizing the friction on 
lateral boundaries are needed to describe F . Finally, the closed equations have been compared to existent continous 
models. Their differences have been discussed: a new term whose form is given should be included to take into account 
lateral boundaries effects. The precise study of the dependence of the three parameters 7, /i and 77 with the beads 
diameters, the beads coefficient of restitution, the beads surface roughness and the cell gap represents interesting 
topic for a future investigation. 
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VIII. APPENDIX 



A. Determination of the volume fraction profile 



3D disordered packing of different volume fractions confined between two plates have been simulated. Numerical 
packing have been constructed by dropping randomly beads in a cell of gap e. Beads are randomly located at a given 
altitude. Their altitudes are then minimized under the two following constrains: the beads should stay in the cell 
and the distance between their center of mass and the center of mass of beads already in the cell should be superior 
to d. Confined packing of volume fraction 0.5 can then be obtained. Confined packing of lower volume fraction are 
then obtained either by contracting beads and cell dimensions without changing the position of beads center of mass 
or by removing randomly a given number of beads. One can then determine the beads actually observed through 
the porthole, i.e. the beads whose center of mass is not hidden by the the presence of the other beads closer to 
the porthole. Both real volume fraction p real and the volume fraction p mew evaluated from the observed beads are 
calculated (see fig. [ll]). p real varies linearly with p mew . For d = 3 mm and e = 7 mm: 

p real = lMp mew (18) 

Direct experimental measurements of p real and p mew just after the filling the drum, i.e. for p real = 0.57 (see section 
III), confirm the value of the ratio p real / p mew . 



B. Determination of the velocity of beads 



The processing to determine beads velocity is the following: each bead is tracked over ten successive frames. As soon 
as the displacement between two successive frames is larger than d, there is an ambiguity and the bead is considered 
as lost. Then, only beads tracked over at least five successive frames are kept. They represent more than 99% of 
all the observed beads. The largest accessible velocity is thus f s d/2 = 1500 mm/s and the error on these velocities 
is about 10 mm/s. The local angle 9 of the surface flow with the horizontal axis is then determined from the mean 
velocity of the flow calculated by averaging the velocity of all beads in all frames. All frames are then rotated by 9 
and divided into horizontal layers one bead diameter wide. The velocity profile is then calculated by averaging the 
velocity of the beads in each layer on 200 frames. The errors made on the velocity profile is then mainly dominated 
by the statistical processing of the data. It is less than 20 mm/s for a 99% confidence interval. 



C. Determination of x, 8 and R 



From experimental H(X) and S(X) profiles, x(X) is calculated via x = J^ P y 1 + (g^(y)) 2 dy where P is the 
upper boundary of the flowing layer (H(Xp) — S(Xp)). The local angle 9(X) is defined as: 9 — atan(dH/dX). The 
procedure to determine R(x) is the following: for a given x, one first finds the corresponding static/flowing boundary 
point Pr of coordinate (X(x), H(x)) in gravity frame. Then one looks for the point Ps of coordinate (X$, S(xs)) at 

the free surface so that PsPh be orthogonal to "e* ' x (numerically, this consists of minimising (X(x) — Xs) cos9(x) + 
{H(x) — S(Xs))sin9(x) with respect to Xs). R(x) is then defined as |PsPh|- 
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• Fig. 1: steady surface flow in a drum at rotation speed O. A thin cascading layer of thickness R (in grey) flows 
on a static bed whose local slope with the horizontal is 9. The transport theorem is applied on the volume 
limited by the vertical slices £(x) and £(x + dx) and resulting equations are written in the curvilinear frame 
(ezjCz). The frame (ex,ez) is the gravity frame located at the center of the drum O. The solid rotation at a 
point M located at the solid/flowing interface is Orri where r is the distance between O and M and ft the unit 
vector perpendicular to olvt 

• Fig. 2: schematic drawing of the experimental set-up 

• Fig. 3: determination of beads velocities. Left: one of the 200 raw frames for = 3 r.p.m.. Frames are binarised 
and beads position is then defined as the center of mass of white areas. Right: corresponding calculated velocity 
field. 

• Fig. 4: left: solid fraction profile (averaged over 1000 frames). Right: velocity profile (averaged over 200 frames) 
at the center of the rotating drum for O = 6 r.p.m.. v x is linear with z in the flowing layer. Inset: v x /v x \ z =o 
vs. the depth showing the exponential creep deformation in the static phase. The straigth line corresponds to 
an exponential fit. In both figures, the errorbars correspond to a 99% confidence interval. 

• Fig. 5: left: velocity profile at the center of the drum for different O: for all these (R, 9) couples, v x is linear 
with z. Right: all these profiles have been translated along the z axis to make their static/flowing interface 
coincide. The velocity gradient 7 is independent of both R and 9. 

• Fig. 6: velocity profile for = 6 r.p.m. at five different locations equally distributed along the free surface. The 
velocity gradient depends weakly on the location in the drum. All these profiles have been tranlated along the 
z axis to make their static/flowing interface coincide. The velocity gradient 7 is consequently independent of 
both d x R and d x 9. 

• Fig. 7: frame processing. Left: one of the 512 raw frames (0 = 3 r.p.m.). Frames are subtracted 2 by 2, binarised 
and gathered by bursts of 25 frames. The resulting frame is smoothed via a dilation-erosion filter. Right: one 
of the 20 resulting frames showing the flowing layer. 

• Fig. 8: typical profiles obtained for O = 3 r.p.m.. Left : profiles H{X) and S(X). Dashed circle correspond 
to rotating drum limits. One can see a large range Cq in the right part of the drum where the layer thickness 
varies while the static/flowing boundary slope almost does not change. Right : profile R(x). 

• Fig. 9: left: mean drift velocity < v x > vs. ^7-R on the range £q. Right: test of the relation v = pR in 
the whole drum. The theoretical profiles Rth (%) and R^ b \x) are calculated from the experimental profile 
H exp (X) using the assumptions (a) v = pR (solid line) or (b) v = pR\J sin 9 (dashed line) (see text for details 
on the procedure). They are then compared to the experimental R exp {x) profile. In all case, p is chosen to 
minimize \R t h ~ R e x P \- 

• Fig. 10: left : variation of F/pjR vs. R for different O. Right: friction force on the lateral plane FB.fr/ P vs. 
R. The solid line corresponds to eq.|l2| (translated vertically of a constant value.): r\ ~ 10 m _1 . 

• Fig. 11: variation of the real volume fraction profile p real vs. the volume fraction profile p mew taking into 
account only the beads effectively seen through the porthole. Data come from numerical simulations. 
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